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ABSTRACT 

We use the new minimum spanning tree (MST) method to look for mass segregation in 
the Taurus association. The method computes the ratio of MST lengths of any chosen 
subset of objects, including the most massive stars and brown dwarfs, to the MST 
lengths of random sets of stars and brown dwarfs in the cluster. This mass segregation 
ratio (Amsr) enables a quantitative measure of the spatial distribution of high-mass 
and low-mass stars, and brown dwarfs to be made in Taurus. 

We find that the most massive stars in Taurus are inversely mass segregated, with 
Amsr = 0-70 ± 0.10 (Amsr = 1 corresponds to no mass segregation), which differs 
from the strong mass segregation signatures found in more dense and massive clusters 
such as Orion. The brown dwarfs in Taurus are not mass segregated, although we 
find evidence that some low-mass stars are, with an Amsr = 1.25 ± 0.15. Finally, we 
compare our results to previous measures of the spatial distribution of stars and brown 
dwarfs in Taurus, and briefly discuss their implications. 
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1 INTRODUCTION 

The Taurus association is a nearby young clust er (at 140 pc, 
with an age of ~ 1 Myr; iKenvon et alj 1994 ). still in the 
process of forming stars from its natal molecular cloud. It 
contains relatively few stars (< 400), of wh ich most are 
contained within several main aggregates (e.g. lGomez et all 
ll993l : lKenvon et al.ll2008h . Star formation in Taurus appears 
to be occurring along three parallel filaments, with the cen- 
tral f ilament coincident on the m ain region of aggregates 
(e.g. lUngerechts fc Thaddeusll 19871) . 

T aurus has a spatial extent of ~ 30 pc (Pall a fe Stahlerj 
2002), and has a low number density compared with, for 
example, the Orion Nebula Cluster. Due to its sparse envi- 
ronment and young age, it is t hought that very little dy nami- 
cal evolution has taken place (|Kroupa fc Bouvierll2003l ). and 
that the observed stars are dir ect signatures of the star for- 
mation process in this region (|Luhmanll2006l ). For this rea- 
son, attempts have been made to quantify the spatial distri- 
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bution of stars and brown dwarfs in Taurus to test various 
formation hypotheses, including those that postulate a dif- 
ferent formation scenario for brown dwarfs ove r stars (e.g. 
iReipurth fc Clarkdl200ll ; iThies fc Kroupall2007l ). 



In this pap er we use the new minimum spanning tree 
(MST) method (| Allison et al.l l2009) to look for differences in 
the distribution of low- and high-mass objects in Taurus. We 
describe the observational sample used in Section [2] before 
presenting the results in Section We compare the MST 
method to other measures of spatial distribution in Taurus 
in Section [4] we discuss our results in Section [5] and we 
present our conclusions in Section [5] 



This is the first in a series of papers in which we will 
discuss the formation of stars and brown dwarfs in Taurus 
by considering the process from prestellar cores to the sub- 
sequent effects of dynamical evolution on the cluster popu- 
lation. 
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Figure 1. A map of the Taurus cluster showing the 361 objects in our dataset. The 20 least massive cluster members are shown by the 
(blue) crosses and the 20 most massive cluster members are shown by the large (red) dots. The a reas of Taurus that are observationally 
complete (surveys by Briccno ct al. 2002; Luhman 2006; Guicu ct al. 2006; Luhman et al. 2010, and references therein) are inside the 
solid lines. 



2 THE OBSERVATIONAL SAMPLE 

Our primary database for the following analysis is a cata- 
logue of 442 Taurus source s compiled by the XEST collab- 
oration (|Giidel et alj|2007l ) as an "input catalogue" for the 
XEST project. This input catalogue was compiled by cross- 
identifying objects between various previous catalogues of 
Taurus members (in part i cular from [Kenvon fc Hartmannl 
1 19951 : iBriceno et all 120021 ; IPalla fc Stahlerl 120021 ) and gen- 
eral all-sky catalogues relevant for pre-main sequence stars. 
Ancillary information such as photometric and spectroscopic 
data, coordinates, masses and ages, was then extracted from 
the individual catalogues, although for some all-sky survey 
catalogues we confined the search for counterparts to within 
a radius of 8 degrees of the position RA(2000.0) = 4 h 25 m 
dec(2000.0) = 25 deg (note that this constraint is irrelevant 
for the identification of Taurus members which relies on pre- 
vious, dedicated Taurus catalogues). Information from SIM- 
BAD and the 2MASS catalogues (essentially spectral types, 
coordinates and photometry) was confined to the areas cov- 
ered by the XEST X-ray exposures (again, this does not 
affect the membership identification relevant for our study). 
A condensed version of the input catalogue for the areas 



covere d by the XEST survey was published in iGiidel et al.l 
|2007l ) where the relevant catalogue bibliography is also de- 
scribed. 

Of these 442 catalogue sources, 293 have a mass esti- 
mate derived from bolometri c luminosit y , Lboi , an d effec- 
tive temperature, T c g, using ISiess et all (|2000l ) isochrones 
with a relative u ncertainty of order of 20 per cent (see 
IGiidel et all 120071 1 . Of the remaining 149 objects without 
a mass listed in the catalogue, 55 have a known spectral 
type . We used this sp ectral type to derive a mass estimate 
from Sicss et al l (2000) isochrones assuming an age of 2 Myr. 
This yielded a total of 328 Taurus members following the 
removal of duplicates. Where available, binary companions 
were included in the sample. Of these 328 objects, 20 do not 
ap pear in the m o re rec ent compilation of Taurus members 
bv lKenvon et alj (|2008l ) and we therefore rejected them from 
the analysis. We will discuss the possible effects of hidden 
binaries and rogue non-members on our results in Section [3] 

The XEST catalogue misses most of the recently dis- 
covered very low mass stars and substellar members of the 
Taurus cloud. We therefore completed the XEST sample 
with the low mass end of the Taurus population taken from 
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iKenvon et alj 's (|2008T ) compilation that lists 382 Taurus 
members. The latter database includes 85 very low mass 
Taurus members not included in the XES T database. Of 
these, only 53 have a spectral type listed in iLuhman et all 
(2010). We used these s pectra l types to derive mass esti- 
mates from ISiess et alj 's l|2000h 2 Myr isochrone. 

Adding these more recent detections to the XEST 
source list eventually yields a catalogue of 361 Taurus mem- 
bers with a mass estimate. We conservatively estimate the 
relative error on the mass to be of order 30 per cent. Alter- 
natively, as a check to the ro bustn ess of our results below, 
we also used lLuhman et alj 's ()2010l ) list of 324 Taurus mem- 
bers with know n spectral t y pes, fo r which we derived a mass 
estimate using ISiess et all 's (|200Ch 2 Myr isochrone. 

We show a map of the Taurus cluster made with our 
data in Fig.[T] The (blue) crosses show the twenty least mas- 
sive objects (all are brown dwarfs) in the cluster, whereas 
the large (red) points show the twenty most massive stars in 
the cluster. Extensive surv e ys of vari o us areas of Taurus by 
iBriceno et al.1 d 19981. |2002|); Luhmanl d200ph: ILuhman et al.l 



1 2003t ); lLuhmanl l|2004 12009 ) ; lGuieu et al.1 (|2006D are shown 
by the black outlines. It is thought that these areas are more 
or less observationally complete , whereas the regions out- 
side of these lines may not be |Luhman et al. 2009, 120101 ; 
iMonin et aill2010l ). 



3 RESULTS 

In this Section we describe the minimum spanning tree 
(MST) method used to quantify mass segregation in clus- 
ters before applying it to sets of objects of similar mass in 
Taurus. 



3.1 The minimum spanning tree method 

Following I Allison et al] f2009), we adopt the minimum span- 
ning tree (MST) method to quantify the level of mass segre- 
gation in Taurus. The MST of a set of points is the path 
connecting all the points via the shortest po ssible path- 
length but which contains no c losed loops (e.g. |Primlll957l ; 

ICartwright fc Whitwor"thll2004l'l. 

We use the algorithm of lPriml (|l957T ) to construct MSTs 
in our dataset. We first make an ordered list of the separa- 
tions between all possible pairs of starfl Stars are then con- 
nected together in 'nodes', starting with the shortest separa- 
tions and proceeding through the list in order of increasing 
separation, forming new nodes if the formation of the node 
does not result in a closed loop. 



3.2 Quantifying mass segregation 

Observationally, 'mass segregation' is a term used to de- 
scribe the central concentration of massive stars in a star 
cluster (the prime example probably being the Trapezium of 
massive stars at the centre of the Orion Nebula Cluster). In 
addition, mass segregation is often used in dynamics to refer 



1 From this point onwards, when referring in general to 'stars' in 
the cluster, we mean 'stars and brown dwarfs', as we are including 
all the objects in the observational sample. 



to the central concentration of massive stars, and the wider 
distribution of low-mass stars caused by energy equiparti- 
tion due to two-body relaxation. 

In this paper we will define 'mass segregation' in terms 
of the relative spatial distributions of stars in a particular 
mass range with respect to other stars in a cluster. This also 
allows us to define 'inverse mass segregation' as an under- 
concentration of a particular stellar mass range with respect 
to the other cluster members. Note that we can apply this 
definition to low-mass stars/brown dwarfs, and by describing 
a population of low-mass stars as 'inversely mass segregated' 
we do not mean that the high-mass stars are necessarily mass 
segregated. 

We find the MST of the ./Vmst stars in the chosen subset 
and compare this to the MST of sets of -/Vmst random stars 
in the cluster. If the length of the MST of the chosen subset is 
shorter than the average length of the MSTs for the random 
stars then the subset has a more concentrated distribution 
and is said to be mass segregated. Conversely, if the MST 
length of the chosen subset is longer than the average MST 
length, then the subset has a less concentrated distribution, 
and is said to be inversely mass segregated. Alternatively, if 
the MST length of the chosen subset is equal to the random 
MST length, we can conclude that no mass segregation is 
present. 

By taking the ratio of the average random MST length 
to the subset MST length, a quantitative measure of the de- 
gree of mass segregation (normal or inverse) can be obtained. 
We first determine the subset MST length, Z S ubset- We then 
determine the average length of sets of Amst random stars 
each time; (Z a verage)- There is a dispersion associated with 
the average length of random MSTs, which is roughly Gaus- 
sian and can be quantified as the standard deviation of the 
lengths (Average) ± ""average • However, we conservatively esti- 
mate the lower (upper) uncertainty as the MST length which 
lies 1/6 (5/6) of the way through an ordered list of all the 
random lengths (corresponding to a 66 per cent deviation 
from the median value, (^average))- This determination pre- 
vents a single outlying object from heavily influencing the 
uncertainty. We can now define the 'mass segregation ra- 
tio' (Amsr) as the ratio between the average random MST 
pathlength and that of a chosen subset, or mass range of 
objects: 



(1 



average 



T 5/6 



T l/6/^subset 



(1) 



A Amsr of ~ 1 shows that the stars in the chosen subset are 
distributed in the same way as all the other stars, whereas 
Amsr > 1 indicates mass segregation and Amsr < 1 indi- 
cates inverse mass segregation, i.e. the chosen subset is more 
sparsely distribute d than the ot h er sta rs. 

As noted by lAllison et all (120091 ). the MST method 
gives a quantitative measure of mass segregation with an 
associated significance and it does not rely on defining the 
centre of a cluster (somewhat impossible for a substructured 
region like Taurus). It also bypasses the various binning 
methods used in determi ning mass segregation through fit- 
ting a density profile fe.g. lAdams et al ] |200ll ; lLittle"fair et al.l 
120031) or tracing the change in mass functio n with radius 
fe.g. lCouliermis et al ]|2004ISabbi" et al ]|2008l ). We shall now 
apply the MST method to look for mass segregation in the 
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Figure 2. The evolution of the mass segregation ratio, AmsRi 
with respect to the TVmst most massive stars in Taurus. Error 
bars show the 1/6 and 5/6 percentile values from the median, 
as described in the text. The dashed line indicates Amsr = 1> 
i.e. no mass segregation. We also show the lowest mass within 
A^mst stars on the top axis. 
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Figure 3. The evolution of the mass segregation ratio, Amsr, 
with respect to the Nmst least massive stars (brown dwarfs) in 
Taurus. Error bars show the 1/6 and 5/6 percentile values from 
the median, as described in the text. The dashed line indicates 
Amsr = 1) i- e - no mass segregation. We show the highest mass 
within Nmst stars on the top axis. 



high- and low-mass stellar (and substellar) populations in 
Taurus. 



3.3 High-mass cluster members 

In Fig. [T] we show the location of the twenty most massive 
objects in the cluster (m > 1.2 Mq) by the large (red) points. 
Several are within the central aggregates, but others are 
located in both the northern and southern Gomez groups 
l|Gomez et al] IT993). In Fig. [2] we show Amsr as a function 
of the number of stars in an MST for the highest mass stars. 
Amsr = 1, indicating no difference between the distribution 
of these stars and other stars, is shown by the dashed line. 

Inspection of Fig. [2] shows that the highest mass stars 
in the cluster are spread more widely than other stars, ie. 
they are inversely mass segregated, with a trough at Amsr = 
0.70 ±0.10. 



3.4 Low-mass (brown dwarf) cluster members 

In Fig. [1] we show the location of the twenty least massive 
objects (all of which are brown dwarfs) by the (blue) crosses. 
Most are concentrated in the central aggregates, but there 
are several in the outlying clumps. We show the calculation 
of Amsr as a function of the number of stars in an MST for 
the low-mass objects in Fig. [3] Again, Amsr = 1, indicating 
no mass segregation, is shown by the dashed line. 

Fig. [3] shows that the distribution of brown dwarfs in 
the cluster is roughly uniform, fluctuating around Amsr = 1, 
with no clear trend towards either mass segregation, or 
inverse mass segregation. There are hints that the low- 
intermediate mass stars may be mass segregated (see Sec- 
tion [33J, an d the brown dwarfs have Amsr < 1, but overall 
the plot is consistent with there being no difference between 
the distribution of low-mass objects and other objects. 



3.5 MSTs for all cluster members 

In a new variation of the MST method, we calculate the 
MSTs for stars as a function of mass. This is achieved by 
taking the MST of a subset of the Nmst lowest-mass ob- 
jects (we take the average of Nmst = 40 objects, rather 
than Nmst ~ 20, to reduce the uncertainties), and then 
sliding through the mass range in steps of 10 objects. For 
example, the first subset contains the 40 lowest-mass ob- 
jects, the second contains the 10 - 50 lowest mass objects, 
the third contains the 20 - 60 lowest mass objects and so 
or0. We then calculate Amsr as before, and plot it in Fig. f3] 
as a function of the highest mass object in each subset. 

It should be noted that in this method the data points 
are not independent of one another, with each data point 
including some of the same information as those in the two 
bins either side. However, if we move through the dataset in 
steps of 30 objects, without any overlap, and compare the 
MST of each subset to random MSTs of 30 objects, the main 
results still hold. 

In Fig. [4] we see again that, when compared to the 
MST of random subsets of objects, the brown dwarfs have 
a mass segregation ratio consistent with unity. Stars with 
masses in the range 0.1 — 0.25 Mq appear to be slightly more 
concentrated (mass segregated), as do stars in the range 
0.45 — 0.8 Mq. In both mass regimes, the mass segregation 
ratio is Amsr = 1-25 ± 0.15. 

Interestingly, stars with masses centred on 0.3 Mq ap- 
pear to have a wider distribution (slightly inversely mass 
segregated), with a trough at Amsr = 0.80 ± 0.10. In Fig. [5] 
we show the location of stars with mass in the range 0.25 - 
0.35 Mq by the plus signs. Most of the stars in our sample 
with this mass have spectral types in the range M2 - M6, 

2 Because we have 361 objects in our sample, the final subset 
contains 41 stars, rather than 40. 
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Figure 4. The evolution of the mass segregation ratio, AmsRi 
as a function of mass for subsets of 40 stars. We plot the highest 
mass object in each subset. The dashed line indicates Amsr = lj 
i.e. no mass segregation. 
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Figure 5. A map of the Taurus cluster, with the positions of stars 
of mass 0.25 - 0.35 Mq shown by t he plus s i gns. T h e deep fields 
survey ed by iBriceno et al. ( 2002); iLuhmanl feoOGl) ; iGuieu et al] 
361). and references therein, are inside the solid lines. 
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in the regime in which the observations may be incomplete 
outside the clu mpy regions of the cluster (|Guieu et al.ll200f3 ; 
iLuhmanllioO^ . 

Our result implies that if there is a deficiency of M2 
- M6 objects, then those that are missing should be lo- 
cated within the clumps, assuming that the anomalous Amsr 
around 0.3 M is a real feature, and that these objects do 
not form via a different mechanism to e.g. objects of mass 
0.2 and 0.5 M . 

Above a mass of ~ 0.9 M , the stars in each subset 
are inversely mass segregated with respect to random stars 
in the cluster, confirming the results shown in Fig. [2] The 
level of inverse mass segregation reaches a minimum value 



of Amsr = 0.70 ±0.10. Whilst this can be said to be a rather 
modest level of inverse mass segregation, it is markedly dif- 
ferent to the MSR for stars with masses of ~ 0.5 Mq. 

3.6 Potential uncertainties 

In this section we briefly discuss the caveats associated with 
our results, namely the main observational uncertainties 
that would affect the resultant Amsr values. 

3.6.1 Mass determination 

The mass determinations for most objects in our observa- 
tional sample are likely to be uncertain by up to 30 per 
cent. It is not possible to directly quantify this in the deter- 
mination of Amsr, as this value is obtained by calculating 
pathlengths between objects, and is not weighted by the 
object's masfl In order to estimate the effect of the mass 
uncertainty on our result, we randomly added or subtracted 
up to 30 per cent of the mass from each object, and then per- 
formed our analysis on this data. From multiple realisations 
of this experiment, we find no significant difference to the 
main result that the most massive stars are inversely mass 
segregated, and the low-mass stars are slightly mass segre- 
gated. However, the inverse mass segregation of objects at 
0.3 Mq is largely erased each time, due to the addition of 
random noise to the mass of each star. The effect of this 
process is place the stars that show strong segregation into 
different bins, diluting the result. 

3.6.2 Binary companions 

We include objects that were lis ted as binary systems in the 
catalogue of IGiidel et al] l|2007h in our analysis. In order to 
test for the effects of close, or hidden binaries that may be 
missing from our data, we performed two experiments on the 
data. Firstly, if an object was multiple, we removed it, and 
its companion(s) from the dataset altogether. This does not 
alter the the results in any way. Secondly, we summed the 
masses of the components and added these to the primary, 
thereby accounting for (and probably overestimating) the 
effects of hidden companions on the mass. Again, negligible 
differences to the main results were found. 



3.6.3 Rogue cluster members 

By comparin g the XEST c a talogu e of IGiidel et al] (|2007t) 
with that in iKenvon et al] ([2008) , we are confident that 
there are no non-members masquerading in our dataset. 
However, should there be any rogue members in our sam- 
ple, they will affect the analysis in 2-D only; i.e. back- 
ground/foreground field stars will not cause an MST length 
to be overly long in the third dimension. Field stars in the 
diffuse regions (outside of the black outline in Fig. [TJ could 
adversely affect the results, but we suggest that the chances 
of this are minimal for two reasons. Firstly, the MST results 
are identical whether we include or exclude the 20 members 

3 An advantage of the MST method is that it does not require 
an absolute mass determination. One can use e.g. absolute mag- 
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of our sample not found in the catalogue of iKenvon et al.l 
J2008I'). Second l y, usi ng the largely independent sample from 
iLuhman etUl |2010h . we also find very similar MST results 
(see Section [4]). This suggests that our observational sam- 
ple would have to change drastically (and that there would 
have to be a significant number of rogue stars distributed 
differently to the cluster members) before the MST results 
are adversely compromised. 

3.6.4 Missing B-type stars 

The fnitial Mass Function (IMF) in Taurus has been the 
subject of much debate. Initially, it was thought that 
Taur us was deficient in both brown dwarfs (iBriceno et al.l 
12002ft and high mass stars (|Walter fe Bovdl Il99lft . This 
contravenes the universality of the IMF, which appears 
the same in most star form ing regions jKroupal 2002; 
iBastian. Covey fc Meyer] 20ld). Re c ently, the discove ry of 
many brown dwarfs 1 Luhmanll2004l ; [Guieu et al.ll2006ft has 
removed the deficit in the low-mass regime. 

However, if one extrapolates the IMF to the high-mass 
regime, there could be up to 4 B-type stars "m i ssing" from 
Taurus l|Walter fc Bovdl Il99ll ). IWalter fc Bovdl |l99lft pro- 
posed that 21 stars in the Cas-Tau OB association were 
related to Taurus. However, 10 of these candidate members 
lie outside the field of view in Fig. [T] and presumably have 
low-mass stars associated with them, for which we have no 
information. To determine the effect of these stars on our 
results, we first added all 21 candidates to our object list, 
before running the MST on this, and a list containing only 
the 11 stars that lie within our field of view. In both cases, 
the net result is the B-type stars are even more inversely 
mass segregated than solar-mass stars. 

In short, if there are missing B-type stars from our ob- 
servational sample, we would expect them to simply rein- 
force our main results. However, we note that in some cases, 
sampling an IMF to populate a low number cluster such as 
Taurus could in principle lead to a deficie ncy in a particular 
mass of object (|Parker fc Goodwinll2007ft . 



3.6.5 Incompleteness in the low-mass regime 

In Fig. [T] the fields for w hich the observation s are thought 
to be entirely complete (jLuhman et al.l l2010l : iMonin et al.l 
|2010| . and references therein) are indicated by the solid lines. 
Outside these regions, it is possible that surveys of Taurus 
may have missed objects, particularly low-mass stars and 
brown dwarfs. Such missing objects may impact upon the 
results of our MST technique. To qualify the potential ef- 
fects of missing objects, we have run the MST on the central 
region only (encompassed by the solid line in Fig. [TJ. The 
results are shown in Fig. [6] The most massive objects have 
a mass segregation ratio Amsr = 0. 81^0 05, which is not as 
extreme a trough as the Amsr = 0.70 ± 0.10 found for the 
whole association. However, in Fig. [1] one can clearly see 
that many of the most massive stars in the association are 
located outside of the central region. If the sparsely popu- 
lated regions in between the central region and the groups 
are more or less complete, then omitting the outlying re- 
gions from the analysis is potentially adding a bias to the 
results because we are no longer considering the entire star 
forming region. 
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Figure 6. As Fig. [2] but for the central region marked by the 
black outline in Fig. [T] The dashed line indicates Amsr = li 
i.e. no mass segregation. We also show the lowest mass within 
^mst stars on the top axis. 



Interestingly, iKirk fc Myers! l|2010ft recently studied the 
sub-groups of stars within Taurus and found that the 
most massive stars in the groups are mass segregated. 
IKirk fc Myersl <|2010h determined the centre of each sub- 
group, and then calculated the offset from the centre for 
each star. They find that the most massive star in each group 
has an offset which is significantly lower than the median. 
We also reproduce this res ult if we calculate Amsr on, for 
example the L1551 group l|Gomez et all 1 1993ft enclosed by 
the black outline at the bottom of Fig. [T] For each of the 
subgroups we find that the most massive stars are mass seg- 
regated. However, the calculated values for Amsr are not as 
robust as those for the whole association, due to low-number 
statistics. We prefer to consider the entire association in our 
analysis, as the sub-groups are interlinked via gas filaments, 
so star formation must be happening on a global scale. 

Additionally, a great deal of observational effort has 
gone into improving the completeness of the entire asso- 
ciation an d it may be that t he data are more or less com- 
plete Ce.g. lGuieu et al.ll200d ; lLuhman et al.ll2010ft . Further- 
more, star formation is recognised to be more prominent in 
filaments, whic h have strong CO signatures and high dust 
extinction (e.g. IPalla fc Stahleri|2002l; ISchmalzl et al llioiol . 
and references therei n). Indeed, Luhman et al. I |2009ft show 
that dust extinction (|Dobashi et alj|2005ft is strongly corre- 
lated with known members, with little evidence of filaments 
elsewhere, suggesting that the stellar census is complete. 

Finally, we note that any theory of star formation in 
Taurus must reproduce the inverse mass segregation we ob- 
serve over a large scale, and the localised mass segregation 
observed bv lKirk fc Myersl (|2010h . We will examine this in 
detail in a forthcoming paper. 

3. 6. 6 Extinction 

A further, related issue to completeness is the variation of 
extinction across the cluster. In the most clustered regions, 
the faintest objects may not be detected due to their being 



On the mass segregation of stars and brown dwarfs in Taurus 7 




Figu re 7. A s Fig. [4] but computed with data provided in 
ILuhman et~aH (2010). The mass segregation ratio, AmsRi is plot- 
ted as a function of the most massive object in each subset of 40 
stars. The dashed line indicates Amsr = li i- e - no mass segrega- 
tion. 



embedded in the gas. As a check, we discarded all objects 
with an extinction A v > 4 and repeated the analysis. Again, 
we find no discernible difference to the results. 



4 COMPARISON WITH OTHER METHODS 
AND DATASETS 

In this section we compare the results of the MST analysis 
of our Taurus dataset with other datasets and with other 
methods that have previously been used to analyse the spa- 
tial distribution of objects in Taurus. 



4.1 Comparison with other data 

In recent work, ILuhman et all (|2010h provided a list of 324 
members of Taurus for which spectral types could be as- 
signed to each object. From these spectral types, m asses 
were inferred using the isochrones of lSiess et alj (|2000| ). As 
an independent test of our method, we repe at the step MST 
analysis in Section 13.51 for the objects in ILuhman et alj 's 
sample and our results are shown in Fig. [7] It should be 
noted that the subsets of objects lie in slightly different lo- 
cations to those calculated using our dataset in Fig. [4] due to 
the fact that there are 37 fewer members overall, and objects 
with similar spectral types are assigned the same masses, 
causing the 'pile-up' of mass segregation ratios at some mass 
values. However, in general, the results are very similar to 
those using our data; the brown dwarfs have Amsr ~ 1, 
whereas stars with masses less than 1 Mq appear mass segre- 
gated, with th e anomalous feature s till prevalent at 0.3 Mq. 
The data from ILuhman et al.1 (|2010t ) are also consistent with 
Amsr = 0.7 (within the uncertainties) for the most massive 
objects in Taurus. 



4.2 The 7Z SS ratio of substellar-stellar objects 

Previous studies into the spatial distribution of brown 
dwarfs in Taurus measured the ratio of brown dwarfs to stars 
for both the whole cluster and the separate aggregates: 



AT(0.02 < m/M < 0.08) 
AT(0.08 < m/M < 10) ' 



(2) 



This ratio has been calculated for the whol e Taurus as- 
socia tion l|Briceno et al.l |2002| ; ILuhmar] |2004| ; iGuieu et al.l 
2006), resulting in a range of values dep e nding on the cho- 
sen dataset. For examp le. iBriceno et al.1 (|2002] ) find TZ BB = 
0.13 ± 0.04 ILuhmanl <|2004l ) finds 1Z aa = 0.18 ± 0.04 and 
IGuieu et al.1 (|2006l ) find TZ BB = 0.23 ± 0.04. IGuieu et al.1 
l|2006h also applied the 1Z BB ratio to the various aggregates 
and concluded that the brown dwarfs are less abundant (by 
a factor of ~ 2) compared to stars in the aggregates than 
for the overall cluster. 

An overall cluster value o f 7^. ss = 0.23±0.04 is consistent 
with the Trapezium cluster (|Briceno et al.1 [20021 ). whereas 
lower values suggest a deficiency in the substellar IMF. How- 
ever. iLuhmanl (feoOot ) argues that the 1Z SS ratio is strongly 
biased by the assignment of spectral type to a particular ob- 
ject (as this changes both the numerator and denominator 
of Eqn.[2j). 

A further, related problem lies in determining the com- 
pleteness of the substellar population. For example, if we 
have 30 brown dwarfs and 220 stars, 1Z BB = 0.14. If a further 
10 brown dwarfs are added to the sample, the ratio of sub- 
stellar to stellar objects becomes TZ aB = 0.18. In other words, 
a normal IMF can appear abnormal simply due to observa- 
tional incompleteness. Such a change to the sample would 
not drastically affect the results of the MST technique, un- 
less the majority of the missing brown dwarfs were spatially 
distributed in a very different fashion to other objects of 
similar mass in the sample. 



4.3 Nearest neighbour distances 

In order to m inimise the perc eived biases associated with 
the 1Z BB ratio, ILuhmar] (|2006t ) adopted the nearest neigh- 
bour distance as a method of quantifying the spat ial distri- 
bution of brown dwarfs in Taurus. In his analysis, ILuhmar] 
(2006) classified objects with spectral type > M6 as brown 
dwarfs, and objects ^ M6 as stellar objects. To ac count for 
the po tential incompleteness in the range M2 - M6 ILuhmar] 
(2006) also made a sub-classification of stars as ^ M2, and 
compared objects with > M6 to bot h ^ M 6 and ^ M2. 

For each object class, ILuhmar] l|2006l ) determined the 
distance to the nearest neighbour. He examined the distance 
from each > M6 (brown dwarf) and ^ M2 (star) to the 
nearest ^ M2; and the distance from each > M6 (brown 
dwarf) and ^ M6 (star - second definition) to the nearest 
^ M6 - see his fig. 14. We repeat his analysis for the dataset 
used here and our results are shown i n Fig. [8] 

We agree with the conclusion of ILuhmanl (2006); the 
distances between brown dwarfs and stars, and stars and 
stars, do not differ by much in our dataset. However, there 
do appear to be subtle variations in the spatial distribution 
as a function of the mass of the object in Taurus (recall 
Figs. [4] and [7]l . These differences are not apparent in the 
nearest neighbour analysis. In Fig. [HI the distributions of 
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Figure 8. The distances to nearest neighbours of stars and brown dwarfs. In (a) we show a distribution of the distances to the nearest 
^ M2 star from: (i) a > M6 brown dwarf (the open histogram with error bars on the left of each bin); and (ii) a ^ M2 star (the hashed 
histogram with error bars on the right of each bin). Each histogram is normalised to the total number of > M6 or ^ M2 objects. In (b) 
we show a distribution of the distances to the nearest ^ M6 star from: (i) a > M6 brown dwarf (the open histogram with error bars on 
the left of each bin); and (ii) a ^ M6 star (the hashed histogram with error bars on the right of each bin). Each histogram is normalised 
to the total number of > M6 or ^ M6 objects. 



nearest neighbour dis tances between any chosen groups of 
objects are identical. iGuieu et alj (|2006l ) find a similar re- 
sult, and both authors found the distribution of stellar and 
substellar nearest neighbour distances to be consistent. 



We therefore caution against using the mean near- 
est neighbour distance to define the spatial distribution of 
brown dwarfs compared to stars in a cluster. If the mass 
function of Taurus is normal, we would expect there to be 
4-5 times as many stars as brown dwarfs in the cluster 
ijAndersen et al.ll200St ). If we calculate the average nearest 
neighbour distance between the brown dwarfs in our sample, 
we obtain a value of 33 arcminutes, compared to a value of 
11 arcminutes between stars. However, this technique is bi- 
ased towards obtaining smaller nearest neighbour distances 
for stars because there are more of these objects in the clus- 
ter than brown dwarfs. Therefore, the stars are more likely 
to be closer to other stars than the brown dwarfs are to 
brown dwarfs. 



If we compare the MST length between brown dwarfs 
to the MST length of random sets of stars, we obtain a 
(largely) unbiased determination of the spatial distribution 
of these objects, and we are also able to pick out the subtle 
differences in the distribution of intermediate mass stars and 
the highest mass stars (see Figs. 2] and 0). 



Finally, we note that other comparisons between the 
MST technique and nearest neighbour distance also find the 
MST to be a more robust de termination of spatial distribu- 
tion l|Gutermuth et al.ll2009h . 



5 DISCUSSION 



We have calculated Amsr (|Allison et al.„ 2009) for stellar and 
substellar objects across the entire Taurus association. We 
find that the most massive stars in the cluster (m > 1.2 Mq) 
are slightly inversely mass segregated with respect to ran- 
dom stars, with a trough at Amsr = 0.70 ± 0.10 (Amsr. = 1 
indicates no mass segregation). This result is unusual in that 
Orion (often considered to be a 'typical' star cluster) dis- 
plays mass segregation of the most massive cluster mem- 
bers (independent of the method used to define mass seg- 
regation), with litt le or no mass segregation below 5 Mq 
(Allison et al. 2009). Currently, the only other cluster to 
have been analysed using the MST method is Trumpler 14, 
and this cluster is similar to Orion in that it displays promi- 
nent mass segreg ation of the most massive stars (> 10 Mq, 
ISana et al.ll2010f ). 

If the data are complete, they suggest that brown 
dwarfs are distributed in a slightly different way to 
most low-mass stars, although within the uncertainties 
the two distributions are fairly similar. However, if brown 
dwar fs form via a different mechanism to low-mass stars 
(e.g. iReipurth fc Clarke! l200ll ; iThies fc Kroupal I2007T ) then 
the observed difference may be real (however se e e.g. 
I Padoan fc Nordlrmdl |2002|. |2004|; IWhitworth et al l 
IStamatellos et al.ll2007l ; iBatei |2009| ; IWhitworth et all 



2007 



2010, 



for arguments that their formation is similar to that of low- 
mass hydrogen burning stars). 

Taking the results of this study at face value lead to 
several conclusions. 

• Firstly, the highest-mass stars in Taurus (m > 1.2 Mq) 
are more widely distributed than average. 
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• Secondly, that brown dwarfs and very low-mass stars 
(m < 0.15 Mq) are distributed randomly in the cluster and 
are not found preferentially either within or outside clumps. 

• Thirdly, that intermediate-mass stars (0.15 < m/Mg < 
0.7) are more concentrated than a random selection of stars. 

• Finally, stars of ~ 0.3 Mq are an exception to the 
concentration of intermediate-mass stars, seemingly signif- 
icantly more widely distributed than stars of even slightly 
higher or lower masses. 

A visual inspection of Figs.[T]and[5]does suggest that the 
first three conclusions are at least plausible, especially that 
the most massive stars are more sparsely distributed. How- 
ever, the finding that stars of ~ O.3M are more sparsely 
distributed than stars slightly more or less massive (Amsr = 
0.8 compared to 1.25) is rather odd and we will return to 
this later. 

Taurus is dynamically young and relatively un- 
evolved. The stella r and gas densities are closely related 
II Gomez et al ] 1 19931 : iMonin et alJl20ld ). and stars are still 
for ming with at leas t 20 prestellar cores found in the clus- 
ter IjKirk et alj|200a ). Therefore, at least to some extent, the 
current positions of the stars follow where they formed. That 
higher-mass stars are found preferentially isolated compared 
to intermediate-mass stars suggests that they formed in dif- 
ferent places. This may reflect how cores fragment, or pos- 
sibly how their masses are distributed. It may be that cores 
that are close together fragment more, forming groups of 
intermediate-mass stars whilst more isolated cores tend to 
form fewer, but larger, stars. Alternatively, perhaps each 
core only produces one or two objects, but that lower-mass 
cores cluster more. (It may be argued that these are two 
are equivalent.) We note that the fragmentation scenario 
should also produce the localised mass segregat i on of the 
sub-groups in Taurus, as found bv lKirk fc Mversl (|2010h . 

Brown dwarfs may be distributed differently to all stars 
of whatever mass. The statistical significance of this result 
is too poor to draw any firm conclusions as the total sam- 
ple size in Taurus is rather small. But this may suggest 
that brown dwarfs form as a different population to stars 
in some way (or that very low- mass cores are distributed 
differently). Strong ejections (e.g. iReipurth fc Clarke! [200 ll ) 
would be expected to provide a fairly strong signature of in- 
verse mass segregation (as dynamics would not have enough 
time to erase much of the signature, iGoodwin et al. 20051) 
and so can prob ably be ex c luded as also found bv lLuhmanl 
(2006); see also IJoergens! |2006). That brown dwarfs are 
not found to be associated with higher-mass stars sug- 
gests that disc fragmentation around larger stars is not the 
form ation mechanism behind most brown dwarfs in Tau- 
rus l|Stamatellos et all 120071 ). We note that gentle liber- 
ation from binaries may give a slightly sparser distribu- 
tion of brown dwarfs when compared to low-mass stars 
(IGoodwin fc Whitworthll2007r i. 

It would seem unlikely that stars of 0.3 Mq would form 
or dynamically evolve in a significantly different way to stars 
of mass 0.2 M© or 0.4 Mq. It is far more plausible that this 
effect is due to incompleteness, or errors in the mass determi- 
nations of these objects. Indeed, the spec tral types that are 
missing, M2 - M6, may be incomplete l|Guieu et al.l l200rj| ; 
iLuhmar] [2006) outside the clumpy regions of the cluster. 
However, for this result to be an artifact of incompleteness 



this particular spectral range must be incomplete inside the 
clumps; more M2 - M6 stars away from clumpy regions will 
make the effect more extreme and not less. For this result to 
be due to incompleteness there must be either (a) more <M2 
and >M6 stars in sparser regions to lengthen the MSTs of 
these types and to lengthen the average MSTs, or (b) more 
stars of M2 - M6 within the clumps. 

Finally, we note that if the masses of all objects in Tau- 
rus were subject to non-systematically change by up to 30 
per cent, then the feature at 0.3 Mq may disappear. Further 
work to better constrain the masses of these objects would 
obviously be desirable. 

We will return to a more detailed theoretical analysis 
of these results in a future paper. 



6 CONCLUSIONS 

We have applied the minimum spanning tree (MST) method 
l|Allison et alj [2009) to search for mass segregation (both 
normal and inverse) in the stellar and substellar populations 
of the Taurus association. To this end, we determine the 
MST length of the 20 least massive stars and compare this 
with the MST lengths of random sets of stars. We repeat the 
procedure for the MST length of the 20 most massive stars. 
The level of mass segregation is then quantified via the mass 
segregation ratio (Amsr, where Amsr = 1 corresponds to no 
mass segregation). 

We also apply a new variation of the MST method to 
compare the MST lengths of subsets of 40 objects to 40 ran- 
dom objects, thereby allowing us to trace the evolution of 
Amsr as a function of object mass. This enables the mass 
segregation ratio of intermediate-mass objects to be calcu- 
lated. 

We determine Amsr for the most massive stars 
(m >1.2Mq) in Taurus and find them to be slightly 
inversely mass segregated (Amsr = 0.70 ± 0.10), 
i.e. preferentially located towards the outskirts of the clus- 
ter. This is unusual in that other star clusters show mass 
segregation of th e most massive stars (|Allison et all 120091 ; 
ISana et alj|201oh . although such clusters are more massive, 
and dense, than Taurus. 

We find that the brown dwarfs in Taurus have a mass 
segregation ratio consistent with no mass segregation, al- 
though we find tentative evidence that intermediate mass 
stars (0.15 < m/Mg < 0.7) show slight mass segregation, 
with Amsr = 1.25 ±0.15. 

These results suggest that brown dwarfs are distributed 
randomly in the cluster, whilst intermediate-mass stars are 
generally concentrated in clumpy regions, and higher-mass 
stars are distributed more widely than average. We note that 
the observations of stellar and substellar objects in Taurus 
may be incomplete for spectral types later than M2, and 
further surveys are desirable in order to determine whether 
low-mass stars are distributed differently to brown dwarfs. 
Whilst incompleteness, especially away from the populous 
well-studied regions may effect our conclusions for low-mass 
stars, it is unlikely that any higher-mass stars are missing 
from the surveys of Taurus and so, unless there is a sig- 
nificant population of low-mass stars away from the known 
clumps, this result is robust. 

Our method avoids the need for the sometimes arbi- 
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trary choice of cluster centre necessary in radially-dependent 
searches for mass segregation. It also directly compares the 
path length between objects of similar mass and random 
objects, rather than the nearest neighbour distance between 
stars and brown dwarfs, or the number ratio of brown dwarfs 
to stars in a particular region and we consider it to be a 
more quantitative measure of mass segregation than previ- 
ous techniques. In a follow-up paper, we will use the MST 
method to compare models of prestellar core fragmentation 
with the observational data (Parker et al. in prep) . 
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